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ABSTRACT 

The light profiles of elliptical (E) galaxies are known to display small systematic 
deviations (0.1-0.2 mag.) from the R^^'^ law. In this paper we show that the senses and 
amplitudes of these departures can be naturally accounted for by a simple distribution 
function constructed on the basis of statistical mechanics of violent relaxation. As a 
consequence, detailed light-profiles can be used to infer about the central potentials of E 
galaxies (the only free shape parameter of our model). Eurthermore, the small deviations 
have recently been shown to correlate with luminosity, C. This observation entails 
a slight breaking of the generally assumed structural homology between E galaxies. 
Using our model we parametrize this broken homology by establishing a correlation 
between a suitably normalized central potential and luminosity. The non-homology 
means that a basic assumption in the interpretation of the Eundamental Plane (EP) 
breaks down, and with it the conclusions derived from it. Instead, by assuming that 
Ai I C is independent of luminosity, we derive a relation akin to the EP directly from 
our correlation. This implies that the EP may have a simple stellar dynamical origin. 
We can reproduce the observed Carlberg-Kormendy relation for the central phase-space 
densities (/c) of E galaxies of identical structure, /c oc but non-homology changes 

it to a much weaker dependence (/c oc C~^'^ for constant Ai/C) which would imply that 
the central phase-space densities of ellipticals are comparable to those of spiral galaxies 
(^ 10~^MQpc~"^(km s ^)~"^). Thus dissipationless merging is consistent with the EP, 
although HST observations, notably the presence of nuclear embedded disks, indicate 
that the assumptions behind our model (isotropy and pressure support) break down for 
the nuclei of faint E galaxies in which dissipative processes seem to be important. 

Subject headings: celestial mechanics, stellar dynamics — galaxies: elliptical and lentic- 
ular, cD — galaxies: fundamental parameters — galaxies: kinematics and dynamics — 
galaxies: photometry — galaxies: structure 



"^Present address: Institute of Astronomy, Madingley Road, Cambridge CBS OHA, UK; jens@mail.ast.cam.ac.uk 



1 



1. Introduction 

Normal, bright elliptical (E) galaxies follow the em- 
pirical R^l^ law (de Vaucouleurs 1948), 



^{R) = + 8.3268 



(1) 



which is a linear relation between the surface bright- 
ness, /J,, and the fourth root of the projected radius, 
R. The effective radius, Rg, is defined as the radius 
of the isophote containing half the total light, and 

= iJ,(Re) is the corresponding surface brightness. 

In itself the universality of this surface-brightness 
profile is a remarkable property, especially in view of 
the very different sizes, fine structures, and forma- 
tion histories of E galaxies. But these systems follow 
additional tight relations. 

The Fundamental Plane (FP) (Djorgovski & Davis 
1987; Dressier et al. 1987) confines the manifold of 
E galaxies to a two-dimensional sheet in the three- 
dimensional space of characteristic luminosity, radius, 
and velocity dispersion. This makes E galaxies a two- 
parameter family. The version of the FP to be used 
in this paper takes the form 



CxRfa'^^: A f« 5 f« 0.8 ± 0.1 



(2) 



(e.g., Schaeffer et al. 1993). This relation is believed 
to be a consequence of the virial theorem, the homol- 
ogy between elliptical galaxies (the R^^^ law), and a 
weak dependence of the mass-to-light ratio, A4/C, on 

M/CxC"; l/a^6±2 (3) 

(Kormendy & Djorgovski 1989; Djorgovski & San- 
tiago 1993). This could be due to dark matter, 
varying IMF or metallicity, etc. (Renzini & Ciotti 
1993; Guzman, Lucey & Bower 1993), but no con- 
vincing theoretical explanation or interpretation has 
been given. For example, many explanations need 
fine-tuning and are expected to have difficulties with 
merger processes which tend to mess up nice correla- 
tions. An intriguing observation in this respect is that 
the FP seems to persist at IR wavelengths where val- 
ues of a 0.1-0.15 are found (Djorgovski & Santiago 
1993; Recillas-Cruz et al. 1991). 

The situation is slightly more satisfactory with re- 
spect to understanding the R^^^ law. Although the 
physics of the violent relaxation process (Lynden-Bell 
1967), which is believed to be an important ingredi- 
ent in E galaxy formation, is still largely unresolved. 



the concept can be used in conjunction with statisti- 
cal mechanics to give a reasonable explanation of the 
i?!/-* law (Hjorth & Madsen 1991, 1993, and refer- 
ences therein). 

Ellipticals are known to exhibit small systematic 
deviations from the R^^^ law, correlated with lumi- 
nosity (Michard 1985; Schombert 1986). For exam- 
ple, Caon, Capaccioli & D'Onofrio (1993) have shown 
that fitting E galaxies with generalized "iji/™ laws" 
leads to a correlation between m and total luminosity. 
This means that the structural homology generally 
assumed for the class of E galaxies is slightly broken. 
Likewise, the end-products of numerical simulations 
frequently show characteristic residuals from the i?^/'* 
law (Londrillo, Messina & Stiavelli 1991). 

Recently Burkert (1993) has treated CCD data of a 
large sample of elliptical galaxies (Bender, Dobereiner 
& Mollenhoff 1988; Franx, Illingworth & Heckman 
1989; Peletier et al. 1990) in a systematic way. The 
combination of a sample of high-quality data with a 
consistent data analysis allows new insight into the is- 
sues mentioned above. The present paper is devoted 
to such an investigation. Our aim is to show that 
the observed departures from the i?^/'* law can be 
understood theoretically essentially as the signature 
of the isothermal sphere. Furthermore, we shall in- 
vestigate to what extent the detected non-homology 
between E galaxies affects the usual interpretation of 
the FP. In fact, we shall argue that the FP may arise 
as a result of this slight breaking of the structural 
homology between E galaxies that nevertheless keep 
approximate i?^/'* law profiles. Finally, in order to 
discuss how mergers may affect the FP we show that 
our approach allows a new assessment of the central 
phase-space densities in E galaxies to the extent that 
these can be considered isothermal and pressure sup- 
ported (e.g., do not harbor embedded disks). 

The paper is built up as follows: In § 2 we briefiy 
review the approach and results of Burkert (1993) and 
introduce a model of elliptical galaxies based on the 
theory of violent relaxation (Hjorth & Madsen 1991). 
The message of this section is that the model is shown 
to be capable of reproducing the observed senses and 
amplitudes of the departures of the i?^/'* law, thereby 
giving them a possible physical explanation. Given 
that this explanation entails a slightly varying struc- 
ture as a function of total luminosity we also describe 
two applications within this framework. In § 3 we 
show that the Fundamental Plane can be understood 
as a result of these slightly dissimilar structures of 
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E galaxies. Finally we compute the central phase- 
space densities of elliptical galaxies in § 4 and show 
that the luminosity-dependent structure of E galax- 
ies has interesting implications for the hypothesis that 
E galaxies may form from dissipationless merging of 
spiral galaxies. In § 5 we discuss the evidence for a 
possible dichotomy between bright and faint galaxies 
in the light of findings of this paper. Conclusions are 
provided in § 6. 

2. Small Deviations from the R^l^ Law 
2.1. Burkert's Analysis 

We first briefiy summarize the results of Burkert 
(1993) (for details, the reader is referred to the orig- 
inal paper). The radius coordinate is defined by tak- 
ing the ellipticity of the observed galaxies properly 
into account. An effective "effective radius", 
determined from the light profile in a self-consistent 
way: From a first estimate of (in arcsec) the best- 
fitting straight line is fitted to the data in the range 
0.6 < X < 1.1, where x is the radius (in units of x^) 
to the power of 1/4. The lower limit in x ensures 
that seeing effects are minimized and the upper limit 
is imposed to avoid fiat fielding errors, sky subtrac- 
tion errors, etc. A new guess of x^ is then determined 
from (8.3268/&)'*a;e, where h is the best-fitting slope 
(see equation (1)). With this updated value oix^ the 
procedure is repeated until convergence (typically 5 
times). 

63 bright E galaxies were studied in this way. 70 % 
of these light profiles showed a characteristic dip at 
X 0.8 being fainter by typically 0.1 mag than the 
best-fitting R^l^ law (for representative observed pro- 
files we refer to Burkert's (1993) paper). The remain- 
ing galaxies typically had a maximum at x 0.9. The 
sense and strength of these departures were quantified 
in the parameter 
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8.3268 



a;G[a;cut;l^l] 



dx 



Xcnt being defined by the minimum (maximum) of 
the bump at x^nt = 0.8 (0.9). The characteristic min- 
imum found in most galaxies gives a negative value of 
bh. In this way Burkert found that luminous galax- 
ies have preferentially negative bh and faint galaxies 
positive bh (see also Fig. 3 below). Thus there is a 
systematic trend in the deviations from the R^l^ law 
that is quantified in a parameter that can be com- 



puted consistently for various models. 

2.2. Model 

In this paper we choose to compare the observed 
profiles with those computed from a theoretical model 
constructed to account for the surface-brightness pro- 
files of spherical galaxies. 

Based on statistical mechanics and a simple vio- 
lent relaxation scenario we have proposed a model for 
the phase-space density in elliptical galaxies (Hjorth 
& Madsen 1991, 1993). It involves a simple two-step 
scenario: First violent relaxation redistributes stars 
in phase space so efficiently that the original predic- 
tion by Lynden-Bell (1967) is applicable. However, 
this only takes place in a finite volume. The second 
step, therefore, is an escape of those particles that 
have gained positive energies during the first violent 
relaxation phase, and redistribution of weakly bound 
particles to larger average radii. 

This scenario leads to a Maxwell-Boltzmann dis- 
tribution in the central parts of the galaxy, / cx 
exp(— E'/cr^); in the outer parts the differential en- 
ergy distribution, i.e., the number of particles with 
a given energy per unit mass, E, (Binney 1982) is 
predicted to be finite, N{E) ^ No, for E ^ 0" but 
zero above the escape energy (for details see Hjorth 
& Madsen (1991)). For numerical confirmations of 
these hypotheses see e.g. Londrillo et al. (1991). 

We have previously discussed the properties of 
the model in terms of its one free shape parame- 
ter, namely the dimensionless central potential = 
— $o/c^- Here $o is the physical central potential 
and (T is a model parameter which in the interest- 
ing range of is approximately equal to the one- 
dimensional central velocity dispersion which is again 
roughly equal to the line-of-sight velocity dispersion 
derived from observations. In these units the model 
predicts 

/T\ ( exp(^) in the inner regions; 
1^ W m the outer regions, 

where n = 4.1 — 4.4, the boundary being approxi- 
mately at the effective radius. In the outer Keplerian 
regions, where $(r) ~ —r~^ this corresponds to a fall- 
off in density as p ~ r~" . Recently, each of these pre- 
dictions have been shown to be in perfect agreement 
with independent numerical #-body simulations of 
merger remnants (Lima-Neto 1994). 

The model has been successful in retrieving the 
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Fig. 1. — The ^ relation. The dashed curve is 
the relation *o = 6.60i/'^''^ used in § 2.2. 

R^l^ law for any ^ 8.5, has given a success- 
ful fit to the small departures from the R^l^ law 
in the light curve of NGC 3379 (Hjorth & Madsen 
1991; Sodemann & Thomsen 1994), and has been 
used to understand radial-orbit instability thresholds 
in dissipationless-collapse simulations (Hjorth 1994). 
In this Section we show that it can also reproduce 
the observed deviations from the R^l^ law found by 
Burkert. 

In the discussion that follows we shall make use 
of a more model-independent dimensionless central 
potential. 



4> 



rh 



$0 



rh 



(6) 



GM GM 

where is the half-mass radius. For example, tp is 
readily computed in numerical A^-body experiments 
(Londrillo et al. 1991). The unique relation between 
^0 and tp is shown in Figure 1. 

This paper mainly deals with the central proper- 
ties of the model (i.e., out to one or two effective 
radii). According to equation (5) this is essentially 
the isothermal sphere. This assertion is confirmed by 
computing the analogous curve to Figure 1 for the 
King model. The latter model also has an isothermal 
core, but a different cut-off in f(E). The correspond- 
ing curve shows very much the same trend. Thus, 
although our discussion will be based on a specific, 
albeit physically well-motivated, model, many of our 
results do not hinge on this choice and our conclusion 
may be generally valid for most realistic models with 
a non-singular isothermal core. 

Some useful expressions for the central properties 



derived from the model (to be used in § 3 and § 4) 
follow below. 

The relation between the half-mass radius and the 
effective radius is rh/Re ~ 1-34 ± 0.01 as for any 
spherical system approximately obeying the R^^^ law 
(Young 1976). Using *o = QMip'^^^ (cf. Fig. 1) valid 
for our model in the relevant range of values of ip (see 
§ 2.3) we obtain from equation (6) the expression 



^2^0.113^^2/^; 

Re 



1-6^1/^^4, (7) 



for the central one-dimensional velocity dispersion, 
which is correct to 1 % for 1.85 < tp < 3.5 and correct 
to 5% for 1.6 < V < 1-85 and 3.7 < V < 4.0. 

A useful expression for the central density is found 
to be 

Pe« 0.0277— 10* (±18%); 1.6 ^ t/^ ^ 4, (8) 

and the central phase-space density, as defined by the 
model, is found to be well-approximated by 



/e « 0.064pe<T-^ (±1%), 



(9) 



(consistent with an isothermal core, fc = (2^)'^^^ pca~^), 
where central velocity dispersions and densities are 
calculated using equations (7) and (8). 

2.3. Results 

It is evident from Burkert's (1993) data that the 
iteratively determined Xg does not give the "correct" 
half-light radius. Repeating the iterative determina- 
tion of Xg for our models as a function of ip confirms 
that Xg may be up to a factor of 2 wrong. However, in 
the present context it is more important that the pro- 
cedure is well-defined: Burkert's observational results 
can be directly compared with model predictions. 
Typical residuals from the R^^^ law (/Umodei — I^r^/*) 
(using Xg, not Rg as the effective radius) for various 
values of ip are shown in Figure 2. 

The lower limit ^ 8.5 mentioned in § 2.2, aris- 
ing from the requirement that the model should give 
reasonable R^^'^ laws, leads to a lower limit on i/" ^ 2. 
Less good but useful i?^/'* law approximations are ob- 
tained down to i/i 1.6 but for these values the iter- 
ative determination of Xg does not converge to a rea- 
sonable value. For these profiles only one redetermi- 
nation of the slope was done. As a consequence, high 
positive values of 6h (^ 0.3) are uncertain. For such 
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Fig. 2. — Computed departures from the R^l^ law 
for i) = 1.71 (dotted curve), V = 2.59,3.07,3.52,4.02 
(solid curves, from top to bottom at x = 1.2), and 
ip = 4.47 (dashed curve). 

small values of i/" a maximum around x 0.85 ± 0.1 
is recovered leading to a positive 6h. 

For slightly larger values of ip 6h approaches 0. 
Then, for a large range of values of ip (from 2.5 to 
4) the observed characteristic minimum at x ^ 0.8 is 
recovered (cf. § 2.1) leading to a negative 6h. The am- 
plitude of the dip also corresponds to those observed 
in E galaxies. Thus, the model describes the small 
deviations found by Burkert satisfactorily. 

For t/" ^ 4 (^0 ^11) this dip becomes less pro- 
nounced and 6h again approaches 0. Furthermore, the 
structure of the model changes dramatically. From 
Figure 1 it is evident that and thus the central 
density of the model (cf. equation (5)) increases dras- 
tically above ip ^ 4. A similar behavior is seen in, 
e.g., the King model (see also Fig. 4-10 in Binney & 
Tremaine (1987)). This is probably an upper limit re- 
lated to the observed upper cut-off in luminosity and 
velocity dispersion for E galaxies. 

This limit is furthermore consistent with numeri- 
cal simulations: Londrillo et al. (1991) found in their 
cold dissipationless collapse simulations the deepest 
central potential to be i/" = 3.20, but dissipation can 
increase the density compared to this value somewhat 
(the simulation leading to the shallowest potential 
had = 1.45 but gave a poor fit to the R^^^ law). 
For the exact ij^/-* law ip = 3.54 (Young 1976). In 
this connection it is interesting to note that from the 
published residuals from the i?^/'* law of Londrillo et 
al. (1991) it is evident that systems with 2 ^ ip ^ 2.4 



have 6h K whereas 6h is negative for 2.4 ^ tp ^ 3.2, 
consistent with our analysis. 

To summarize, for 2.5 < ip < 4 the model repro- 
duces the decrease in brightness at x ^ 0.8 relative 
to the R^^^ law observed in 70% of bright E galaxies 
(6h < 0). For 1.6 < i/" < 2 it can also reproduce the 
observed increased brightness relative to the R^^^ law 
(6h > 0) for the remaining galaxies. 

3. The Fundamental Plane 

The question now is whether ip can be determined 
from photometry of elliptical galaxies. From the 
model we can calculate 6h as a function of ip and 
Burkert found a dependence of 6h on C To exam- 
ine whether there is a relation between ip and C we 
now incorporate the FP in the discussion. 

For the purposes of our discussion, the FP relation 
(2) can be written in the simplified form 

£(X (i?e,72)T; 0.7^7^0.9, (10) 

without any decisive loss of generality. Combining the 
observational FP relation (10) with equation (7) we 
obtain 



Two simple interpretations of this equation are 
possible: 

(i) If Ip is independent of C then A4/C (x C" . 

(ii) If Mj C is independent of C then ip"^!^ cx C" cx 
M". 

Normally structural homology (i) is assumed, i.e., 
that all galaxies have similar structures. In this case 
equation (3) is recovered. This is a reasonable hy- 
pothesis given that most E galaxies are well-described 
by an exact R^l^ law. 

The second possibility is to take the systematic de- 
partures from the R^l^ law as an indication of dissimi- 
lar structures, corresponding to a ip-C correlation. In 
the simplest case (ii) Ai/C is independent of £. This 
hypothesis can be tested by comparing the resulting 
value of a with that determined observationally. 

The structure-luminosity trend can be visualized 
in a bh-Ms diagram. Burkert's (1993) relation (his 
Fig. 8) is not unique, but this is not surprising given 
that the observations are difficult and that much of 
the fine structure of elliptical galaxies (e.g., disks, 
shells, dust, environment etc.) will undoubtedly spoil 
such a relation somewhat. What is interesting is that 
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Fig. 3. — Circles are bh values for a subsample of 
galaxies (Burkert 1993) with C4 < and e < 0.3 con- 
taining no dust as a function of total B magnitude, 
Mb ■ The thin vertical line is the approximate bound- 
ary between galaxies with m > 4 (left) and m < 4 
(right) found by Caon et al. (1993). The curves are 
relations obtained using Mb = — 2.51ogi/'^''(^") — K 
with (a = 0.09, A' = 12.5) (dotted curve), (a = 
0.12, K = 15.5) (solid curve), (a = 0.18, A' = 18.5) 
(dashed curve). All values of Mb are scaled to 
h = 0.5. 

there does seem to be a trend: For Mb ^ —19.5 (Mb 
is the total absolute blue magnitude) 6h is positive 
whereas 6h < for all galaxies with Mb ^ —22. Fur- 
thermore, there seems to be a lower boundary on 6h. 
These are the trends that we can hope to reproduce. 

In order to minimize the potential effects arising 
from possible contributions from disks in the E galax- 
ies a subset of Burkert 's sample is plotted in Figure 3. 
These data points only include galaxies with "boxi- 
ness parameter" C4 < (e.g., Peletier et al. 1990) and 
ellipticity e < 0.3. Furthermore, only galaxies show- 
ing no signs of dust or other peculiarities are plotted. 
It is seen that there is still a clear trend although 
there are much fewer data points (22). 

In a related study Caon et al. (1993) found that 
fitting i?^/™ laws in general leads to m > 4 for 
Mb ^ —19.3 and m < 4 for fainter galaxies. Sim- 
ilarly, Young & Currie (1994) have found that dwarf 
elliptical galaxies all have m < 4 with a conspicuous 
correlation between m and Mb extending the findings 
of Caon et al. (1993). Although its physical meaning 
is less clear, the parameter m is similar to 6h in that 
higher m leads to higher ip (Ciotti 1991) and that 



m — 4 and 6h have opposite signs. The findings of 
Caon et al. (1993) and Young & Currie (1994), which 
are based on a larger sample of galaxies covering a 
larger range in absolute luminosity, thus strengthen 
the conclusions derived from Figure 3. 

Assuming that the possible range of ip for bright E 
galaxies suitably described by the model is 1.6 ^ i/" ^ 
4, the value of the exponent a (cf. equation (11)) is 
quite well-determined if the allowed range in ip must 
span the observed magnitude range. In Figure 3 we 
show a good fit to the 6h vs. Mb data obtained by 
letting a = 0.12t° °^ (7 = 0.88 ± 0.03) under the 
hypothesis (ii) of constant Ai/Cg. This value of a 
is in agreement with observations, especially in the 
infrared (~bolometric) where Djorgovski & Santiago 
(1993) found a ^ 0.13 and Recillas-Cruz et al. (1991) 
found a 0.15 for galaxies in Virgo and a 0.09 
in Coma. Observations in infrared light are prefer- 
able since they more directly probe the old relaxed 
stellar population and have the advantage that line- 
blanketing effects are negligible. The observed larger 
values of a at shorter wavelengths is to a great deal 
probably due to a metallicity effect (Renzini & Ciotti 
1993). 

Assuming a correlation between the structure of 
galaxies (parametrized by ip) and their luminosity we 
have thus shown that the model leads to a FP-like 
relation consistent with observations even with a con- 
stant mass-to-light ratio. 

4. Central Phase-Space Densities 

To look into the possible origin of this FP rela- 
tion, in particular the effects of merging, we next 
compute the central phase-space densities of E galax- 
ies. Dissipationless merging cannot increase the maxi- 
mum (central) phase-space density (Tremaine, Henon 
& Lynden-Bell 1986). Therefore, accurate determina- 
tion of central phase-space densities may have impor- 
tant consequences for formation theories of elliptical 
galaxies, in particular the merger hypothesis (Vedel 
& Sommer-Larsen 1990). 

Unfortunately, determination of central phase-space 
densities of elliptical galaxies is difficult. Previous de- 
terminations have used the core-fitting method (Rich- 
stone & Tremaine 1986) which is based on a fit of an 
isothermal model to the central parts of the galaxy. 
However, ground-based observations are strongly af- 
fected by seeing and most galaxies are known to have 
non-isothermal cores, with a significant excess of light 
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near the center (Kormendy 1985a; M0ller, Stiavelli & 
Zeilinger 1993). For example, recent HST (Hubble 
Space Telescope) data of the cores of ellipticals in- 
dicate very high central values, presumably due to 
central black holes (Crane et al. 1993). 

We here propose a different method based on the 
indirect determination of the central potential using 
the light-profiles in the regions un-affected by seeing 
and the central cusp. In the adiabatic model the 
phase-space density is unaltered by a black hole so 
our assumption of an isothermal core may be ade- 
quate for the present purposes. 

To compute physical phase-space densities (in units 
of MQpc~^(kms )~^) using the expression for fc 
given in equation (9) we shall resort to an approxi- 
mate relation between Rg (in parsec) and Mb ■ For 
purposes of illustration we choose to combine the 
Hamabe-Kormendy (1977) relation i^gB = 2.941ogi?e-|- 
11.93 (valid for a distance modulus for the Virgo clus- 
ter of (m — M) = 31.3 and a Hubble constant h = 0.7; 
Hq = 100/ikms~^ Mpc~^, cf. Capaccioli, Caon & 
D'Onofrio (1992)) with Ms = -5 log i?e-F//e_B -24.05 
which follows from the i?^/'* law (Young 1976), to 
get Mb = -2.061ogi?e - 12.12 2.94 log(0.7//i), i.e.. 
Re oc Cg"^^^. Equations (7) and (8) then lead to 

/e OC lO^r'^B'''- (12) 

This relation clearly illustrates that there is a dra- 
matic effect of galaxies having ip dependent on lumi- 
nosity: For the usual assumption that all ellipticals 
have a fixed structure (ip constant) one recovers (up to 
a prefactor depending of the choice of ip) the relation 
of Carlberg (1986) who found fc oc using data 

of Kormendy (1985b). However, using our Mb — ip 
relation with a = 0.12 a much weaker dependence 
of luminosity, fc oc is obtained, and the values 

are found to be quite low, ^ 10~®MQpc~^(km s ^)~^. 
These effects are illustrated in Figure 4. 

Before discussing this result a word of caution is 
appropriate. The results shown in Figure 4 were de- 
duced assuming a scaling relation between Rg and 
Mb ■ The existence of the FP and the distribution of 
galaxies on the FP means that such a relation is only 
approximate and contains a lot of real scatter. What 
our discussion does illustrate, however, are the poten- 
tially important effects that the varying structure of 
ellipticals may have. Furthermore, assuming different 
relations between Rg and Cb than the one adopted 
here results in different fc^MB relations (cf. Fig. 4) 
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Fig. 4. — Central phase-space density log (in units 
of M0pc-3(kms"^)-3) versus Mb ■ Filled circles are 
exact values for the model using Rg oc jCb"^^^ (see 
text). The solid curve is calculated using the approx- 
imations given in the paper. The dashed line is the 
approximate relation log 0.6Mb + 4.5. The dot- 
ted line is the much steeper Carlberg-Kormendy re- 
lation (Carlberg 1986) valid if structural homology 
is assumed. The dashed-dotted curve is similar to 
the solid curve but has been derived using Rg oc Cy^^ 
(Guzman et al. 1993) and 5 - = 0.95. /Cg = 10 
and h = 0.5 have been assumed in this plot. 

but the qualitative effects of the varying structure are 
unaltered. 

5. Discussion 

Taken at face value, the results that we find in- 
dicate that most ellipticals may have central phase- 
space densities that are comparable to those of spi- 
rals for which fg K (lO-5-lO-^)M0pc-3(kms"^)-3, 
cf. Carlberg (1986) and Lake (1989). Notice in par- 
ticular that the central phase-space densities of faint 
elliptical galaxies are predicted to be about two orders 
of magnitude smaller than predicted by the Carlberg- 
Kormendy relation. The reason why the small devia- 
tions, i.e., the central potential, has such a dramatic 
infiuence on the measurement of can be traced back 
to the fact that any model with a sufficiently deep cen- 
tral potential can lead to the R^l^ law. Assuming a 
unique profile for all galaxies corresponds to chosing 
a specific central potential. As we have shown here, 
this may be erroneous and can lead to substantially 
different derived results for the central properties of 
the galaxies. 
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As an illustration of this effect consider the data 
sample of Kormendy (1985b), upon which the Carlberg- 
Kormendy relation is based. The central projected 
surface brightnesses (and thus the central phase-space 
densities) of the galaxies were obtained by perform- 
ing a seeing deconvolution of the observed light pro- 
files. Following Schweizer (1981) it was assumed that 
the cores of the galaxies were described by a King 
model with c = 2.25 (corresponding to ip ^ 9.5, or 
Mb = —22.5 for a = 0.12, i.e., a quite deep central 
potential). But as M0ller et al. (1993) have demon- 
strated the assumptions underlying the seeing decon- 
volution are very model dependent and may give mis- 
leading results. 

It can be seen from Figure 4 that the Carlberg- 
Kormendy relation and the one found in this paper 
(equation (12)) agree quite well for the bright galax- 
ies. These galaxies have the smallest seeing correc- 
tions (Kormendy 1985b) and are well-described by a 
c = 2.25 King model because they have deep dimen- 
sionless central potentials, ip. The large discrepancies 
occur for the faint unresolved galaxies for which the 
seeing corrections (and thus the systematic uncertain- 
ties) are largest. Furthermore, according to the find- 
ings of the present paper, the faint galaxies have cores 
that are better described by King models with con- 
siderably smaller values for the concentration (down 
to c 1.6). The seeing deconvolution would therefore 
overestimate the value of the central surface bright- 
ness and thereby the central phase-space density for 
the faint galaxies. 

Still, the fact that our relation falls significantly 
below the Carlberg-Kormendy relation is somewhat 
disturbing. Although difficult to compute without a 
useful model for the distribution function recent HST 
observations (Crane et al. 1993) indicate that the 
central phase-space densities are indeed much higher 
than predicted by Figure 4, especially for faint galax- 
ies. However, our model is only valid for E galax- 
ies that are pressure supported and isothermal and 
there seems to be a dichotomy between bright and 
faint galaxies, the dividing line being approximately 
as indicated in Figure 3: Bright galaxies in general 
are resolved with the HST (Kormendy et al. 1994), 
do not contain nuclear disks and have boxy isophotes 
(Nieto, Bender & Surma 1991), belong to the 'bright 
group' of Capaccioli et al. (1992), are pressure sup- 
ported (Davies et al. 1983), have 6h < and m > 4, 
and have central phase-space densities less than those 
for spirals. On the other hand, the faint E galax- 



ies in general are unresolved with HST, they often 
contain nuclear dust and gaseous/stellar disks, some- 
times with signs of star formation, belong to the 'ordi- 
nary group', are rotationally supported, have 6h > 
and m < 4, and have phase-space densities larger than 
those of spirals according to the Carlberg-Kormendy 
relation. 

It is therefore possible that the discrepancy be- 
tween the inferred high central phase-space densities 
from HST and our estimates in Figure 4 is a measure 
of the amount of dissipation that has occurred, i.e., 
bright galaxies may have formed essentially dissipa- 
tionlessly through a merging hierarchy whereas (the 
nuclei of) faint E galaxies involved significant dissipa- 
tion. Star formation and material in disks may signif- 
icantly raise the inferred central phase-space density. 
It should be noted, however, that the 'underlying E 
galaxy' may well have a much lower central phase- 
space density and be well-described by our model fur- 
ther out (where 6h is measured). Also, the central 
properties affected by dissipational processes may be 
unimportant for the gross dynamical evolution of an 
E galaxy (excluding the nucleus) due to the very little 
mass fraction involved. In this case the dissipationless 
scenario may still be valid. 

Thus, to the extent that E galaxies can be consid- 
ered pressure supported and isothermal our proposed 
interpretation of the FP is consistent with a hierar- 
chical dissipationless merging scenario for the forma- 
tion of E galaxies in which the mass is increased and 
the central phase-space density is decreased in each 
step. This conclusion is supported by recent indepen- 
dent #-body simulations of Capelato, de Carvalho & 
Carlberg (1995) who find FP-like relations and devia- 
tions from the i?^/'* law quite similar to those found in 
this paper from a merging hierarchy of galaxies with 
constant mass-to-light ratio. 

Correlations between luminosity and other observed 
dynamical quantities may also be envisaged in view 
of the dichotomy (or, perhaps, gradual change with 
luminosity) discussed above. For example, additional 
stellar dynamical effects can be a changing velocity 
structure with luminosity, such as rotation (Prug- 
niel & Simien 1994) or velocity anisotropy (Merritt 
1988). In fact, Burkert (1994) has shown that all 
galaxies in his sample with anisotropic velocity dis- 
tributions, (v/a)* < -0.33, have -0.25 < 6h < 0. 
Assume for example that an increase in radial ve- 
locity anisotropy is correlated with luminosity. For 
an Osipkov-Merritt Jaffe model (Binney & Tremaine 



8 



1987) we find that the line-of-sight velocity dispersion 
within an aperture of i?e/5 or i?e/10 increases by a 
factor of 1.2 when going from an isotropic velocity dis- 
tribution (faint galaxies) to a very radial velocity dis- 
tribution (bright galaxies) with an anisotropy radius 
Tajr-ji = 0.1. Over 5 absolute magnitudes in luminos- 
ity this variation in the observed line-of-sight velocity 
dispersion would translate into a value of 7 less than 
1 (cf. equation (10)) if the galaxies were assumed to 
have identical velocity distributions. For our example 
this corresponds to a = 0.09. 

6. Conclusion 

In this paper we have used a simple model for 
the phase-space density distribution in E galaxies de- 
rived from statistical mechanics of violent relaxation 
(Hjorth & Madsen 1991) to give a possible explana- 
tion of the small deviations from the B}l^ law found 
by Burkert (1993) in a sample of 63 galaxies. In this 
interpretation, the observed departures from the B}l^ 
law are essentially the signatures of the isothermal 
sphere, p cx exp(^), cf. equation (5) and Lima-Neto 
(1994). This view is also consistent with the numer- 
ical simulations of Londrillo et al. (1991). Further 
progress along this line could be obtained by analyz- 
ing a larger sample of high-quality profiles (e.g., Caon 
et al. 1993) in a consistent manner that can be car- 
ried out also for a set of models. To bridge the entire 
interval from giant to dwarf ellipticals it will be neces- 
sary to consider an alternative parameter describing 
the light profiles since the iterative determination of 
bh does not converge for faint galaxies. 

Our model and the observed correlation between 
the deviation of the light profile from a pure B}l^ law 
(as measured by bh, cf. equation (4)) leads to a rela- 
tion between the central potential and the luminosity 
for E galaxies. 

Following this interpretation we have shown that 
the Fundamental Plane (FP) for the manifold of el- 
liptical galaxies can be understood as the effect of a 
slightly changing structure of elliptical galaxies while 
they keep B}l^ law profiles with deviations of only 
0.1 mag (cf. Fig. 2). Adopting a universal mass- 
to-light ratio, Ai/C, for all galaxies we have deter- 
mined the FP "A^/£" exponent (cf. equation (3)) 
a 0.12 from these deviations in agreement with 
independent observational determinations in the in- 
frared. We have thus demonstrated that it is possi- 
ble, using simple stellar dynamics, to get a FP rela- 



tion and that an Ai / C variation with luminosity is 
not needed on purely dynamical grounds. In this pic- 
ture the FP arises as a result of the slightly broken 
homology between E galaxies. Recently, a similar in- 
terpretation has been suggested to be viable for the 
FP correlations found for the Galactic globular clus- 
ter system (Djorgovski 1994). 

To test this view it would be interesting, if possi- 
ble, to establish the FP using IR array photometry. 
Likewise, to test for possible evolutionary effects (dy- 
namical or chemical) it is important to investigate the 
FP correlations of E galaxies at significant redshifts. 

As an alternative to the usual core-fitting tech- 
nique we have determined the central phase-space 
densities of elliptical galaxies using their detailed light 
profiles as an indicator of the depth of the central po- 
tential. This method has the advantage of being see- 
ing independent. We find that the central phase-space 
densities are comparable to or less than those of spi- 
ral galaxies (^ 10~®MQpc~^(km s "'^)~^) and follow a 
relatively weak dependence of luminosity, cx C'j^^ . 
This indicates that ellipticals as well as spirals can 
merge into ellipticals without much dissipation and 
can occupy the FP without having a special IMF 
or distribution of the dark-matter. We attribute the 
high central surface-brightness of E galaxies observed 
with HST to dissipational processes, e.g., the forma- 
tion of embedded nuclear disks or AGN. 

Finally, we point out that our results have been 
reached in several tempi: First, we have made the case 
that the deviations from the R^l^ law and their corre- 
lation with luminosity are evidence for non-homology 
of elliptical galaxies. Second, we have argued that 
the central potentials of ellipticals can be determined 
from our model. Third, in our discussion of the FP 
and phase-space densities, we have assumed a con- 
stant mass-to-light ratio (case (ii)). This assumption 
constitutes the simplest case, but in reality one could 
have a combination of non-homology and a variation 
of mass-to-light ratio with luminosity. 

JH acknowledges discussions with Inger j0rgensen 
and financial support from the Carlsberg Foundation. 
Audi Burkert kindly provided the data points upon 
which Figure 3 is based. This work was supported in 
part by the Danish Natural Science Research Coun- 
cil (SNF) and the Theoretical Astrophysics Center 
(TAG) sponsored by the Danish National Research 
Foundation. 



9 



REFERENCES 

Bender, R., DSbereiner, S., k MSllenhoff, C. 1988, 
A&AS, 74, 385 

Binney, J. 1982, MNRAS, 200, 951 

Binney, J., & Tremaine, S. 1987, Galactic Dynamics 
(Princeton: Princeton Univ. Press) 

Burkert, A. 1993, A&A, 278, 23 

Burkert, A. 1994, MNRAS, submitted 

Caon, N., Capaccioli, M., k D'Onofrio, M. 1993, MN- 
RAS, 265, 1013 

Capaccioli, M., Caon, N., k D'Onofrio, M. 1992, MN- 
RAS, 259, 323 

Capelato, H. V., de Carvalho, R. R., k Carlberg, R. 
G. 1995, ApJ, submitted 

Carlberg, R. G. 1986, ApJ, 310, 593 

Ciotti, L. 1991, A&A, 249, 99 

Crane, P. et al. 1993, AJ, 106, 1371 

Davies, R. L., Efstathiou, G., Fall, S. M., Illingworth, 
G., k Schechter, P. L. 1993, ApJ, 266, 41 

de Vaucouleurs, G. 1948, Ann. Astrophys., 11, 247 

Djorgovski, S. 1994, ApJ, in press 

Djorgovski, S., k Davis, M. 1987, ApJ, 313, 59 

Djorgovski, S., k Santiago, B. X. 1993, in Struc- 
ture, Dynamics and Chemical Evolution of Ellipti- 
cal Galaxies, Proc. ESO/EIPC Workshop, ed. I. J. 
Danziger, W. W. Zeilinger, k K. Kjar (ESQ Conf. 
and Workshop Proc. 45) (Garching: ESO), 59 

Dressier, A., Lynden-Bell, D., Burstein, D., Davis, R., 
Faber, S., Wagner, M., k Terlevich, R. 1987, ApJ, 
313, 42 

Franx, M., Illingworth, G., k Heckman, T. 1989, AJ, 
98, 538 

Guzman, R., Lucey, J. R., k Bower, R. G. 1993, MN- 
RAS, 265, 731 

Hamabe, M., k Kormendy, J. 1987, in Structure and 
Dynamics of Elliptical Galaxies, ed. T. de Zeeuw 
(lAU Symp. 127) (Dordrecht: Reidel), 379 

Hjorth, J. 1994, ApJ, 424, 106 

Hjorth, J., k Madsen, J. 1991, MNRAS, 253, 703 

Hjorth, J., k Madsen, J. 1993, MNRAS, 265, 237 

Kormendy, J. 1985a, ApJ, 292, L9 

Kormendy, J. 1985b, ApJ, 295, 73 



Kormendy, J., k Djorgovski, S. 1989, ARA&A, 27, 
235 

Kormendy, J. et al. 1994, in ESO/OHP Workshop on 
Dwarf Galaxies, ed. Meylan, G. (Garching: ESO), 
in press 

Lake, G. 1989, AJ, 97, 1312 

Lima-Neto, G. B., 1994, A&A, submitted 

Londrillo, P., Messina, A., & Stiavelli, M. 1991, MN- 
RAS, 250, 54 

Lynden-Bell, D. 1967, MNRAS, 167, 101 

Michard, R. 1985, A&AS, 59, 205 

Merritt, D. 1988, AJ, 95, 496 

M0ller P., Stiavelli M., & Zeilinger W. W. 1993, in 
Structure, Dynamics and Chemical Evolution of 
Elliptical Galaxies, Proc. ESO/EIPC Workshop, 
ed. I. J. Danziger, W. W. Zeilinger, & K. Kjar 
(ESO Conf. and Workshop Proc. 45) (Garching: 
ESO), 131 

Nieto, J.-L., Bender, R., & Surma, P. 1991, A&A, 
244, L37 

Peebles, P. J. E. 1993, Principles of Physical Cosmol- 
ogy (Princeton: Princeton Univ. Press) 

Peletier, R. F., Davies, R. L., Illingworth, G., Davies, 
L. E., & Cawson, M. 1990, AJ, 100, 1091 

Prugniel, Ph., & Simien, F. 1994, A&A, 282, LI 

Recillas-Cruz, E., Carrasco, L., Serrano, A. P. G., & 
Cruz-Gonzales, I. 1991, A&A, 249, 312 

Renzini, A., & Ciotti, L. 1993, ApJ, 416, L49 

Richstone, D. O., & Tremaine, S. 1986, AJ, 92, 72 

Schaeffer, R., Maurogordato, S., Cappi, A., & Bernardeau, 
F. 1993, MNRAS, 263, L21 

Schombert, J. M. 1986, ApJS, 60, 603 

Schweizer, F. 1981, AJ, 86, 662 

Sodemann, M., & Thomsen, B. 1994, A&A, in press 

Tremaine, S., Henon, M., & Lynden-Bell, D. 1986, 
MNRAS, 219, 285 

Vedel, H., & Sommer-Larsen, J. 1990, MNRAS, 245, 
637 

Young, P. 1976, AJ, 81, 807 

Young, C. K., & Currie, M. J. 1994, MNRAS, 268, 
Lll 

This 2-column preprint was prepared with the AAS ]J^T^]X 
macros v3.0. 



10 



